Optical Simulation of Yang-Baxter Equation 
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In this paper, several proposals of optically simulating Yang-Baxter equations have been presented. 
Motivated by the recent development of anyon theory, we apply Temperley-Lieb algebra as a bridge 
to recast four-dimentional Yang-Baxter equation into its two-dimensional counterpart. In accordance 
with both representations, we find the corresponding linear-optical simulations, based on the highly 
efficient optical elements. Both the freedom degrees of photon polarization and location are utilized 
as the qubit basis, in which the unitary Yang-Baxter matrices are decomposed into combination of 
^—s ■ actions of basic optical elements. 
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Yang-Baxter equation (YBE) originated in solving the one-dimentional <5-interacting models [l[ and the statistical 
O ■ models on lattices 12, |3| . The importance of YBE was further realized as a starting point for the quantum inverse 
scattering method [4|. It is well-known that YBE plays important role in solving the integrable models in quantum 
field theory and exactly solved models in statistical mechanics ([3| and references therein). In quantum field theory, 
YBE describes the scattering of particles in (1+1) dimensions. The essence of YBE is to factorize the scattering 
of three particles into successive two-body scattering processes. YBE also plays an important role in completely 
integrable statistical models, whose solutions can be found by means of the nested Bcthc Ansatz [5]. 
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Due to its importance, YBE deserves thus to be tested experimentally. Measuring the spectrum of spin chain, 
which is at the root of YBE, one can learn the structure of spinon and thus check the factorization of YBE. For 

£f-\ ' instance, Hcisenberg spin- 1/2 chain model has been probed experimentally through neutron scattering experiments 
and the spectrum coincides with the calculation based on YBE [6]. However, YBE only provides sufficient condition 

t"***" ' for the prediction of spectrum. So these experiments should be viewed as indirect check of YBE. The direct verification 
is still an open question. 
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In order to keep the paper self-contained, we first explain the basic formula of YBE. The Yang-Baxter matrix R 
is a N 2 x iV 2 matrix acting on the tensor product space V ® V, where N is the dimension of V. Such a matrix R 
satisfies the YBE 

u w II -r- 1) w w u 7/ -r- 1) u 

Rx2(u)R2s( —jr 2 )R 12 (v) = R 23 (v)R 12 ( + )R 23 (u). (1) 

1 + p z uv 1 + p^uv 



where Ri 2 = -R <S> 1, R23 = 1 (8 R, u and v are spectral parameters, fi^ 1 — ic (c is light velocity). Take the two 
spin-1/2 particles as an example. Acting on such a system, R is a 2 2 x 2 2 matrix whose matrix elements are R a b,cd, 
a, b, c, d =1 (spin up), J, (spin down). If spin is conserved, some matrix elements may vanish. The physical meaning of 
R(u) is two-particle scattering matrix depending on the relative rapidity tanh _1 (/?M). When we change the spectral 
parameters as 0u = (1 — x)/(l + x), (3v = (1 — y)/(l + y) and j3(u + v)/(l + j3 2 uv) = (1 — xy)/{\ + xy), we obtain 



another ordinary form of YBE 



Ri2{x)R 23 {xy)Ri2{y) = Raa(y)Ri2(xy)Raa(x), (2) 



i.e. the spectral parameter in the middle i?(xy)-matrix being the product of the neighborhoods' spectral parameters. 
The asymptotic limit R(x — ► 0) = b satisfies the braid relation 

612&23612 = b 2 3bi 2 b 2 3. (3) 
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FIG. 1: (a) Diagrammatical interpretation of braid operators The labeled particle worldlines, orienting from bottom 
to ceiling, are unaffected by smooth deformations in which the lines do not intersect. Each crossing means a 

scattering of two particles, including permutation process as the special case, (b) The braid relation J3|) 

&i2^23^i2 = &23&12&23- (c) Skein relation, (d) The unknotted loop. By skein relation we can eliminate all the 

crossings and get a linear combinations of the Kauffman brackets for various disjoint unions of unknotted loops. 

This relation is diagrammatically represented in Fig. {!]). For a given matrix b satisfying ([3]) we can retrieve the 
corresponding ^(a;)-matrix via the procedure of Baxterization or Yang-Baxterization [7|. This procedure depends 
on the numbers of independent eigenvalues of matrix b. In particular, when a braid matrix b has two independent 
eigenvalues Ai and A2, the corresponding ^(a;)-matrix obtained via Yang-Baxterization takes the form 



R(x) = p(x)(b + x\i\2b 



(4) 



where p(x) is normalization factor. For statistical models on lattice, the elements of R(x) should be positive-definite, 
since they are related to the Boltzmann weights. However, as we will see below, there is no such a restriction for the 
application to quantum entangled states. 

In the recent years there is a new development to connect the braid matrix, as well as YBE, with the entangled 
state [3, M, LL1J LLU, LL2, LL3|, LbS, LL5| • The basic idea comes from the Bell states having the maximal entanglement degree. 
For a two-qubit system, Bell states are defined by 



i$±) = _ ( itT)±m)), 

|* ± ) = 4=(|U>±|4T». 
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The Bell states are connected to the natural basis l^o) 
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(We refer the reader to [13J for more details about the short notation in the first line equation.) L. H. Kauffman et 
al have shown that the matrix W is nothing but a braid matrix satisfying (J3]) by recognizing V as two dimensional 
complex vector space to hold a single qubit of information [9(. Further, it was found that W can be extended to 
matrix b such that [13| 
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where real parameter </> is time-dependent flux. Yang-Baxterizing this braid matrix b, one can define a new state with 
arbitrary entanglement degree, as follows 



\*{x,q)) = R(x,q)\* Q ), R(x,q) 
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(8) 



where R(x, g)-matrix satisfies YBE ([2j and meanwhile determines the evolution of the initial state \^q) to \*f?(x,q)) 
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given the time-dependent q = q(t). In terms of the new variable cos© = (1 + x)/(y/2(l + x 2 )), R(x, q) can be recast 
to 



R{Q,<t>) = 
It is interesting to observe that R(Q, </>) satisfies to the relation 



V 



cos 6 







e - ^sin0\ 







COS0 


- sin 9 







sin 9 


cos 8 


e 1 ^ sin 


e 





cos 6 / 



(9) 



R(Q,(f>) = cos 91 + sin QM, M A = -1, (10) 

i.e. as the extension of the Euler formula. When = 7r/4, -R(0, 4>) reduces to 6 in Q, i.e. yielding the maximum of 
entangled states. When takes other values the state |^(0, (f>)) = i?(0, 4>)\^o) processes a continuous entanglement 
degree determined by which is usually less than the maximum [l3j |. Suppose only the flux <f> depends on t and 
under the adiabatic approximation, then we can obtain the Berry phase related to YBE [131 ]. 

The theory sounds reasonable, but why we prefer to choose a Yang-Baxterized matrix R(Q, 4>) as the unitary 
evolution is just an assumption. Especially, one may doubt the necessity of introducing the third particle to describe 
two-particle entanglement. We have to present a practical scheme to test the YBE in the framework of quantum 
information. Fortunately, there have been popular optical operations for the quantum gates [161 ] that are available to 
experimentally test the validity of such a YBE. The motivation of this paper is to propose a linear optical simulation 
of YBE based on the highly efficient optical elements. 

Both the braid matrix B in ([7]) and the .R(x)-matrix in (|9|) act on the tensor product space V (8) V and thus 
have the 4-dimensional (4D) representation. The entangled state in ([5]) requires that the optical simulation of the 
corresponding ^-matrix involves the universal entangled gate, i.e. CNOT gate. In principle CNOT gates make use of 
measurement-induced nonlinearity and are still of low efficiency by means of optics [la . Il7| ] . The situation becomes 
worse when several sequent CNOT gates are involved. We have to find alternative ways to avoid this difficulty. 
Fortunately, as we will see in Sec|TTl the 4D Yang-Baxter matrices have two-dimensional (2D) counterparts which are 
unitary and have much simpler realization by means of linear optics. 

The paper is organized as the following. In SecJITl we first prove the equivalence between 4D braid matrix and 2D 
braid matrix, then find the Yang-Baxterization procedure for 2D braid matrix. Based on this theoretical assertion the 
optical test of 2D YBE will be presented in Sec lIIII A direct test of 4D YBE is shown in Sec lIVI The conclusion is 
made in the Sec[V] The relationship between the basic for 2D braid matrices and 4D ones will be given in Appendix 
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II. TWO TYPES OF YBE'S 

In the topological quantum computation theory, the 2D braid behavior under the exchange of anyons [18| has 
been investigated based on the v ~ 5/2 fractional quantum Hall effect (FQHE) 11511. Motivated by this interesting 
application of braid relation in anyon theory, we will nest Temperley-Lieb algebra [lj| into 4D YBE and reduce it to 
2D YBE. Here we briefly present such an equivalence between these types of YBE's. 

Let us first recall the braid behavior in v = 5/2 FQHE. Quasiparticles in such a system are often called Ising 
anyons or SU(2)2 states, which satisfy non-Abelian fractional statistics. There are three types of anyons, which can 
be called 0, 4, 1. When two anyons become close together while other anyons are much farther away, these two anyons 
can be treated as a single particle whose quantum numbers are obtained by combining the original quantum numbers. 
For SU(2)2 states, such a formation of new anyons obeys the following fusion rules 

ixi=0+l,ixl = i,lxl = 0, 

0x0 = 0, 0xi = i 0x1 = 1. (11) 

(These fusion rules are analogous to the decomposition rules for tensor products of irreducible SU(2) representations, 
but have an important difference that 1 is the maximum spin.) Note that there are two different fusion channels 



for two \ anyons. As a result, when four | anyons fuse together to give 0, there is a two-dimensional space of such 
states. This can be done by dividing the four \ anyons into two pairs. Both pairs either fuse to or to 1 then fuse 
the resulting anyons together to form 0. The orthogonal basis states read |15| 
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In the middle fusion chains (called conformal block), the internal edges are subject to the fusion rules at each trivalent 
vertex. In such conformal block basis, exchanging anyons is identified as braiding in Fig. [TJ From the conformal basis 
to the Kauffman graph in the right-hand sides, Jones- Wcnzl projector operators have been applied, i.e. 
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where d — \J~2~ in present case. By means of skein relation in Fig. [IJ one can introduce the braid operators A and B 
which have nontrivial action on lei) and \e2) 
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Thus their matrix representations in the basis (|ei), \e2}) are given by 
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and they satisfy 2-dimensional braid relation 
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ABA = BAB. 
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We emphasize that (fTB"]) act on the basis (|ei), \e2))- It is worthy of noting that the "crossing" in (Ti"4"l) means the 
usual 4x4 braid matrix, satisfying ([3]). These braid matrices A and B are unitary and have a natural realization by 
linear optics, as we will see in Sec lIIII 

In order to generalize the above procedure from braid relation to YBE, we nest Temperley-Lieb algebra [19j into 
4D YBE and surprisingly reduce it to 2D YBE. Detailed calculations will be given later (see Appendix |A"|) . Briefly, 
acting on the subspace spanned by |ei) and |e 2 ), the 4D YBE (TT|) will reduce to the corresponding 2D YBE 

Mu)B L U , + J )A(v)=B(v)A( U+ J )B(u), 
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where p(u) is normalization factor and e = ±1. Since these matrices A(u) and B{u) are unitary, it is obviously easier 
to optically simulate this 2D equation than the previous 4D edition. For the convenience of experimental check, we 



further introduce the transformation 



1 + P 2 u 2 + 2iej3u 
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then we obtain the following form of SU(2) matrices 
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In terms of this new parameters, the solution of i?(#)-matrix takes the form 
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(20) 



Though it takes the similar form with ([9]), the parameter 9 has different meaning from the parameter 0. One may 
wonder where is the missing parameter cf> in 2D YBE. It actually survives in the basis |ei) and |ea) which 2D YBE 
should act on (see the explicit form of the basis in Appendix |A|) . 
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FIG. 2: Schematic setup for simulating either side of 2D YBE (jTT)) by means of polarization qubit. The angle 
parameters 0i are determined by (|25[) while the relations between 9i and u, u are determined by (|24[) . 



These two unitary matrices can be optically realized with the aid of the " polarization qubit" or " location qubit" 
of a single photon. In the following we present two experimental setups to simulate YBE by means of these two ways. 



III. OPTICAL SIMULATION OF 2D YBE 



The single-photon representation of qubits plays an important role in linear optical computation [2CJ, l2l|, 122J • The 
key is that a single photon can be utilized to act both as polarization qubit and as location qubit. For the former 
one encodes the qubit in the photon's polarization with the corresponding transformations simulated by wave plates, 
such as half- wave plates (HWPs) and quarter- wave plates (QWPs). For the latter one encodes the qubit in the 
single-photon paths with the corresponding transformations simulated by beam splitters (BSs), phase shifters (PSs), 
and mirrors. As a result, universal unitary gates (1-qubit or 2-qubit) can be realized by means of either of these two 
kinds of qubit transformation, or by combining both [20l . l2ll ] . 



The two unitary matrices in (fl9|) can be optically realized with the aid of the "polarization qubit" or "location 
qubit" of a single photon. In the following we present two experimental setups to simulate YBE by means of these 
two ways. 



A. Using Polarization Qubit to Simulate 2D YBE 



We first recall the action of a QWP upon the basis states of the polarization qubit [2 

U Q (S) 
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(21) 



where <7j are Pauli matrices and 5 is the angle between the QWP axis and the vertical direction. Then the action of 
a HWP upon the basis states of the polarization qubit is given by 



Uh(S) = U 2 {8) = -i 



cos(2J) sin(2£) 
sin(25) -cos(25) 



(22) 



As an analogue with Euler rotation, the sandwich configuration of one HWP and two QWPs enables one to perform 
any unitary changes of the photons polarization state [23[ . Particularly, for the case in (|19p , we obtain 



A(0) = U Q (j)U H (-j + -)C/ Q (J), B(0) = U Q (~)U H C-)U Q fy. 



(23) 



By this decomposition it is straightforward to design the experimental setup for simulation of 2-dimensional YBE 
(TTT)) . As Fig[2] shows, a suitable series of QWPs and HWPs with different direction angles in succession will simulate 
the left-hand side (LHS) of YBE while another series will do the right-hand side (RHS). In Fig[2J the relation (Q2 
requires the angle parameters 9i in the LHS satisfy 
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This also holds for the angle parameters appearing in RHS of YBE (|17p . What we should measure in experiment 
is to analysis the actions of transformation of both sides given the input states with the same angle parameters, for 
example, by means of quantum state tomography [24| . 
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FIG. 3: Gates acting on location qubit basis: (a) Beam splitter Ubs- (b) Phase shifters Up S (£,) (left) and Up S (£) 
(right), (c) Hadamard gate H by a 50:50 BS and two — it/2 PS. (d) Mach-Zchnder interferometer Umz as universal 
1-qubit gate. For simplicity, each pair of — 7r/2 phase shifters accompanied with every beam splitter as in (c) is not 

shown in (d). 



B. Using location qubit to simulate 2D YBE 

When we take photon location paths as qubit basis, the unitary transformations can be achieved by means of BSs 
and PSs. We follows the notations in [20| which are different from those in [2lj|, especially the opposite definitions of 
location qubit lead to distinct actions of mirror. First, we list the actions of several elementary gates on the location 
qubit basis (see FigOJ 
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Based on these gates, a Mach-Zehnder interferometer (FigJ3[d)) can realize arbitrary U(2) group element [2(J,[2l|]. Note 
that in Fig[3Jd) each pair of —n/2 phase shifters accompanied with every beam splitter is not shown for simplicity. We 
hold this convention hereafter, so each BS should be taken as Hadamard gate. The unitary action of Mach-Zehnder 
interferometer is given by 



U MZ = Upsi^HU^rrllpsi^Upsi^HUpsifo) 
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(27) 



where the total phase <I> = ^1 + ^2 + ^1+952 and phase difference A = if2 — ^>i- Through Mach-Zehnder interferometer, 
we can perform the action of operators A and B in the (|23|) with angles correspondences as 



A{9) = U M z(<p2 = (f>i = 0, 

B(9) = Umz{^>2 = —ipi = 
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Then we come to the whole optical setup to simulate both sides of 2D YBE, (jTTJ), as shown in FigdJ The angle 
parameters obey the same relation in 
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FIG. 4: Schematic setup for simulating either side of 2D YBE (fl7|) by means of location qubit. (a) Simulation of 
LHS. (b) Simulation of RHS. Each pair of — tt/2 phase shifters accompanied with every beam splitter is not shown. 

The relations of different parameters are refereed to (|2"5|) . 



IV. OPTICAL SIMULATION OF FOUR DIMENSIONAL YBE 



The YBE in four dimension representation fTJ) is equivalent to 

(A(0i) ® i 2 ) ■ (i a ® A(0a)) • (A(e 3 ) ® i 2 ) = (i 2 ® £(0 3 )) • (#(fe) ® 12) • (i 2 ® fl(0i)), 



(29) 



where (9^ satisfy the relation (|25[) and -R(#) takes the form in (j2"0|) . The R(8)-matiix can be decomposed into the 
combination of elementary gates [25j]. The case 9 = is trivial. When 6* = 7r/4 or 37r/4 (we restrict 9 £ (0, 27r)), R{&) 
reduces to the braid matrix ((7|) and thus equivalent to one CNOT gate, as L. H. Kauffman et al first pointed out [3]. 
When 9 takes other values, it can decomposed as follows 



R{6) = (Fi y 2 ) • CNOT2 • (V 3 ® 14) • CNOT2 • (F 5 <g> F 6 ), 
where V* € J7(2) and CNOT2 gate is given by 
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Below we focus on the general case, i.e. 6 =/= 0, 7r/4 or 37r/4. Recently, S. S. Bullock and his coauthors [26[ has 
developed a criterion for determining the number of CNOT (or equivalently CNOT2) gates to simulate a given 
transformation. By this criterion, one can check that the decomposition (|30|) is optimal, i.e. R(9)-m&trix admits a 
quantum circuit using 2 CNOT gates (see Appendix [C|) . 
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FIG. 5: (a) A circuit for simulating left-hand side of 4-dimensional YBE ([5]), with s and p denoting location and 
polarization qubit channels respectively, (b) the decomposition of -R(#)-matrix, see ([30]). 



Detailed calculation (see Appendix |C|) further shows Vi € SU(2) for the present case. More explicitly, we have 
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By this explicit decomposition (|3H)) . the design of the circuit to simulate YBE is straightforward, as shown in Figj5] 
The difficulty lies in realization of the CNOT gates. As is described in previous section, photon can carry either 
"polarization qubit" or "location qubit" . If we only use the former, CNOT gates are possible for photons in principle 
using measurement-induced nonlinearity [271 ] - However, currently they are still low-efficient and experimentally ex- 
pensive [la, Ll3| ■ For the present status of linear optics experiments, it was shown that the success probability of an 
array of n CNOT gates can be made to operate with a probability of p = (^)™ +1 [23|- The above decomposition of 
R(d) takes 2 CNOT2 gates. For each side of YBE (0, we have to deal with 6 CNOT gates at the same time. Thus 



the success probability is p = (-|) 7 ~ 4.57 x 10 4 , which makes the practical simulation extremely difficult. This is 
the reason why we did map the 4D YBE to 2D YBE. 

Photon used as location qubit will help reduce the above difficulty when we do small-scale quantum calculation. 
The key lies in the high efficiency of BSs, PSs and wave plates. In Fig[5]we have to deal with three qubit, then two 
schemes are available: one polarization qubit channel plus two location qubit channels, or all location qubit channels. 
We focus on the former since it save one optical way and use less number of beam splitters. 

In FiglUa), three channels are designed to be location (s), location (s) and polarization (p) qubits from top to 
floor. For polarization qubit, using the QWPs and ITWPs, the unitary matrices Vi can be decomposed into 



Vi = U Q (^)U H ^)U Q C-, 
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For the location qubit, as the previous section shows, a Mach-Zehnder interferometer can simulate these Vi matrices. 
We summarize the results in Fig|6fa) and (b). 

The rest come to two types of CNOT2 gates: one is between location qubit and polarization qubit, the other is 
between two location qubits. For the former, it can be achieved by a polarizing beam splitter where the location qubit 
is flipped or not conditionally on its state of polarization, as shown in Fig(6^c) [2fj]. For the latter, the two location 
qubits in FiglBJd) correspond to the first and second number of the binary representation of the location of a single 
photon, respectively. Thus the corresponding CNOT2 is simulated by simply swapping the labels of path 1 10) and 
|11) (see FigEfd)). 

Gathering all the elementary gates, we finally arrive at the whole scheme for optically simulating the LHS of 4D 
YBE, as shown in FigJTJ In order to get the optical setup to simulate the RHS of 4D YBE, we can apply the formal 
equivalence between two hand-sides in the YBE (JlJ by cycling the indices 1 — ► 2, 2^3, 3 — *• 1 and exchanging the 
parameters u <-> v. The equality will be confirmed by means of tomography given the same input on each sides. 
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FIG. 6: (a) The realization of U\ acting on location qubit with 0i = - 
U2 acting on location qubit with <f>i = %, <fe= ^%: <Pi = 4^ an d ¥>2 



tt+4> 



and 6-2 = 



5TT + 4 



(b) The realization of 

2 = — jr 1- ) fi := ; ^T^ arLcL ^2 = y 4 "" ■ (c) CNOT2 gate using a polarizing 
beam splitter with the polarization and location being the control and target qubit, respectively, (d) CNOT2 gate 
between location qubits, which is achieved by swapping the labels of output paths |01) and |11). 
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CONCLUSION 



We have presented several proposals to optically simulate Yang-Baxter equations. According to the development of 
theoretical analysis, Yang-Baxter equation in two-dimensional representation and in four-dimensional representation 
can be uniformed with the aid of Temperley-Lieb algebra. In both representations, we have found the corresponding 
linear-optical realizations, based on the highly efficient optical elements, i.e. half-wave plates, quarter-wave plates, 
beam splitters, phase shifters, and mirrors. Both the degrees of freedom of photon polarization and location have 
been utilized as the qubit basis. In each kind of basis, the unitary Yang-Baxter matrices have been decomposed 
into combination of actions of basic optical elements. The developed proposals, in principle, are able to be used to 
directly check Yang-Baxter equation. We remark that the test of 2D YBE is, in fact, to provide an optical simulation 
of two-component anyons associated with FQHE. The anyon behavior is now a hot topic [15], but in our knowledge 
there has not been a scheme to test by using an optical simulation. 
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APPENDIX A: REDUCTION OF 4D YBE 

Analogue with the mapping from 4D braid relation to 2D braid relation in SecHIl we here give the calculation of 
their Yang-Baxterized edition. This is equivalent to Yang-Baxterize (p~6|) . First we recall the involved YBE 

Ri2(u)R 23 ( )R 12 (v) = R 23 ( V )R 12 ( —?- )R 23 ( U ). (Al) 

1 + p z uv 1 + p^uv 

What is interesting is that this four-dimensional YBE (|A1[) admits the celebrated Temperley-Lieb algebra (TLA). 
Actually set 

R 12 (u) = ai(u)l6 + 6i(w)/7i 2 , , A2 ^ 

R,23(u) = a 2 (u)t 6 + b 2 (u)U 23 , 

and suppose U satisfying TLA [19| 

U 2 = dU, U 12 U 23 U 12 = J7i2, U 23 U 12 U 23 = U 23 , (A3) 



where d is the loop in FigfTJ taking the value of \/2 in our case. The coefficient functions ai{u) and bj(u) in (|A2[) 
are determined by the associated YBE (|A1|1 . Consider one simple but important case: ai(u) = a 2 (u) = a{u) and 
bi(u) — b 2 (u) — b(u). We easily get 

II -4- 11 II -I- 11 

[a(u)b{v) + b(u)a(v) + d b(u)b(v)]a(-—^-) = [a(v)a{u) - b(v)a(u)]b(f—^-). (A4) 

J. | LL U J- | Hi U 

H) has the solution 



a(u) = p(u), b{u) = p(u)G(u), G{u) = ^ (e = ±1). (A5) 

V 2(1 + p z w z — ziepu) 

On the other hand, in accordance with (fT3)) and (fT4| . the new basis |ei) (i = 1, 2) are introduced from the definition 
of the operator U 

U12 |ei) = C/34 |ei) = d|ei) , 
U l2 |e 2 ) = t/ 34 |ea) = 0, 

t/23 |ei) = U u |ei) = i(|ei) + Vd 2 - 1 |e 2 », (A6) 

a 



Thus we have 



U 23 |ea) = C/14 |ea) = -^ (|ei) + Vd 2 -l\e 2 )) 



Rii{u) |ei) = [oi(«) + d 61 (u)] |ei) , 
Ri2(u) |ea) = ai(u) |e 2 ) , 

&3(«) lei) = [a2(«) + ^M] |ex) + ^^(u) |e 2 ) , ( A7 ) 

^2 _ I d 2 - 1 

R2s(u) \e 2 ) = b 2 {u) |ei) + [a 2 (u) + ——b 2 (u)} |e 2 ) . 
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So |ei) and \e<2) span the _R-invariant subspace. Then it is natural to define the matrix elements A(u)ij = (ei|^i2(u)|ej) 
and S(u)y = {ei\R2z(u)\ej) . 

Combing the above results we obtain the explicit form of A(u) and B(u), 



1 + (3 2 u 2 + 2ie/3u 







A{u) = p{u) ( l + /J2 U 2 _ 2ie(iu 

1 



B(u) = 



p(u) 



1 + (3 2 u 2 2iej3u 
<■-,:■! ■■>;, I-',, \ 2ie(3u 1 + /3 2 u 2 



1 + (3 2 u 2 - 2ie(3 2 u 
and importantly, in accord with the braid relation (flo]) . they satisfy the two-dimensional (2D) YBE 

The corresponding unitary matrix £/ with d = v^, which satisfies TLA in (|A3|) . takes the representation 



V2 



/ 1 ig" 1 \ 
V -iq 1 / 



(A8) 



(A9) 



(A10) 



Here the important factor i before matrix M distinguishes U from braid operator in ([7]). In terms of new parameters 
as in (fT8|) . the explicit form of R(6) takes 



R{B) = a(u) + b(u)U = p{u){t 4 + G{u)U) = 



( cos 9 e - ^ sin 6» ^ 

cos 6» -sin 6» 

sin6» cos6> 

\-e^sin0 cos j 



(All) 



By setting a — p and b — pp, where p — ^ {~d i Vrf 2 — 4) with d = v2, i.e. p = — exp(±i7r/4), we regain A and 
B matrices as in (|T5|) . which satisfy the braid relation (fTB|) . This result can be also obtained through the "light-cone" 
limit of (|A1|) by setting three arguments in ^-matrices to be equal, i.e. 



u + v 



u = v = 



l + ^uv 1 

which is satisfied by either u = V = or u = v = ft -1 . Under the limit in (|A12|) . A(u) and B(u) reduce to ([T5 
We know that by taking Q = iq~ x = i (i.e. Q A = 1), U matrix given by (|A10|) becomes 



(A12) 



U(Q = i) = 



1 



1 i 
1 i 



V2\ -i 1 

-z 1 



(A13) 



which is the transformation matrix for the Bell States [l0( ■ We thus conclude that the 4-dimensional entangling braid 
matrix (|A13|) and the 2-dimensional braid matrix (fTl))) can be uniformed by acting the TLA operator on different 
dimensional basis. The 4-dimensional basis can be (| f|"),| |1), | J.1"), | JJ.)), whereas the logic qubit basis read 
(|ei), |e2)) . Conversely, the latter can be expanded in terms of four-spin states and thus relate with the 4-dimensional 
basis. In order to find this correspondence, we firstly rewrite Uij as a form of projectors 



where 



Uij = v^GV'y'K^M + \<Pij){<Pij\), 



\1>ii) = i(l TT>« +e-^'\ Uh), \tp tj ) = -J=(| U)ij - i\ iT>«), 0' 



37T X 



(A14) 



(A15) 
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It is interesting that both \ipij) and \<Pij) are maximally entangled states for two spins, i.e. Bell states. With the aid 
of (|A6|) . we arrive at 



l ei ) = R , i ,0 (1^12)^34) + ^^12) |V34)), 

|e a ) = - 7 i=((l - ive^ihaMu) ~ *(1 - ^ 1 e^')l^ 3 )|^i)) - Id), (A16) 

Vi + kl 

with v be an arbitrary coefficient. Detailed calculation is shown in Appendix IBl 

Briefly, in the invariant subspace spanned by |ei) and |e2), the R matrices satisfying 4D YBE (|Alj) will reduce to 
2D representation, A(u) and B(u) 1 with the corresponding 2D YBE (fT"7|) . 



APPENDIX B: CALCULATION FOR |e*> 

Here we give the details of calculation for |e») (z = 1, 2). We start from 

U 12 \e 1 )=U 3i \e 1 )=d\e 1 ), (Bl) 

U 12 |e 2 ) = C/34 |e 2 ) - 0, (B2) 

U 23 |ei) = J7i4 |ei) = ~(|ei) + ^d 2 - 1 |e 2 )), (B3) 

a 

U 23 |e 2 ) = C/i4 |e a ) = ^EI (| ei ) + V^I |e 2 )). (B4) 

a 

Generally, |e*) can be expanded into the linear combination of Bell states 

|^) - -^(| Tt)ii±e-^'| ||) y -), |p±) = -L(| U)« T»| it)«)- (B5) 

Due to the project form of t/y 

^•=V2(|^+)(^| + |4.)(^.|), (B6) 

one can get general expression of |e») from (|B1|) and (|B2[) 



|ei) =011^2)^34) + ^1^2)1^34) + a 3^1 2 )l^3 f 4) +^1^12)1^34). 

|e 2 ) =as|V'r2)IV'3"4) + a 6|l/ , r2>l^34> + a 7 \<Pu) I ^34) + a s|^i2> W 3 a) ■ (B7) 

We further notice that (|B3[) and (|B4|) indicate the symmetry of exchanging pair indices 23 <-> 14 for |e<). Taking of 
this symmetry and noticing the minus sign in the expression of ip^ and ip^, we further simplify |ej) into 

|ei) =Ol|^12>l^34> + a 4|Vr2)IV?34)> 

|e 2 ) = a 5 | ^2)^34) +«8|vr 2 )l^34)- ( B8 ) 

Now we set d = v2 in (|B3p and (|B4[) . which means C/ 2 3|ei) = f/i4|e2)- This is the only one condition that further 
determines the coefficients a* (i = 1,4,5,8). What we proceed is detailed expansion 

u 23 \ ei ) = V2{\{\ it) + e -*' I ii))«TT I + e*' (IT I) + i(| U) - i\ lt))((U I + *<IT l)} 23 

{yd TT) + e-«\ |i))i 2 (| TT) + e- # 'l ii)) 3 4 + y (| TI) - *| IT))i 2 (l TI) - i| IT)) 34 } 

= l^ + 3 ){y(l TT) + e-^'| li»i4 + yHI II) - te*'| TOM 
+ l¥&H T (e-*'| TI) + ie-«'\ IT))i4 + y(-| IT) + *| Tl))i4} 

= ^ ~ ie ^' ^Mta^u) + ^(aie-^' +ia4,)\cpi 3 )\ipu), 
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U u \e 2 ) = V2{\{\ IT) + e-*' | ii»((lT I + e**' (ii I) + \{\ Ti) - *| IT))((U I + m l)} 23 

{yd TT) - e-^'| ii))i 2 (| TT> - e"**'| ii)) 34 + f(| U) + i| 4T»i 2 (l Ti) + i| IT>) 34 } 
= h/^Hyd TT) + e-^'| ii))i4 + y^| ii) + ie^| TT»u} 

+ l^3>{y (-e-^'| Ti) - »e-*'| it))i4 + y(-| it) + »| U))u} 
= i(a 5 +ze^a 8 )|V 2 + 3 )|^+) + ^(-aee^' + ia 8 )|^+ )|^ 4 >- (B9) 

Because of the orthogonality of Bell states, we get the relation between coefficients 



a\ - 


- ie 


i<j>' 


(24 = 


a 5 


+ ie 


'■<t>' n 






a\e 


—i(j>' 


+ 


ia,i - 




-a 5 e" 


-»y + j 


ag, 




=> 


a 5 = 




-ie** 


04 


, a 8 


= — ie" 


-20' 


ai 



(BIO) 



Setting ai = — ^^^^= and 04 = — : , we finally arrive at 



v\ 2 



1 



|ea) = - 7 =^=(^^'|Vr 2 )lV'34) +e- i<A W 2 >»)- (BH) 

They are in deed equivalent to (|A16[) . The arbitrary parameter v represents a certain degeneracy between the 
components of |ej) with respect to the actions of Uij. From the process of calculation, we can view D93 and Uu as 
the entanglement swapping operators on the pair-entangled states |e»), in accord with the results in [13| . 



APPENDIX C: DECOMPOSITION OF fl(0) 



Here we give the proof of the decomposition (f30|) based on the work of S. S. Bullock and his coauthors [25|, l26l ] . 
They have developed the following criterion 

Proposition 1. An operator u £ 5/7(4) can be simulated using no CNOT gates and arbitrary one-qubit gates from 
SU(2) iff x[M u )] = (x ± l) 4 . Here 7(14) = "u(er a ® a y )u T (a y (g) cr y ), w T denotes the transpose and x[ff] = det[x/ — g] 
denotes the characteristic polynomial of g. 

Proposition 2. An operator u € 5(7(4) can be simulated using one CNOT gates and arbitrary one-qubit gates 
from SU(2) iff x[l{u)] = (x + i) 2 {x - if. 

Proposition 3. An operator u £ 5/7(4) can be simulated using two CNOT gates and arbitrary one-qubit gates 
from SU(2) iff tr[ 7 (u)] = is real. 

Direct calculation shows that the case 9 = satisfies Proposition 1 while the case 9 = ir/4 or 37r/4 satisfies 
Proposition 2. Since x[t(-^(^))] = (l + x 2 — 2x cos 29) and ti[y(R(9))] = 4 cos 2$, so Proposition 3 confirms that 
R(9) generally admits a quantum circuit using two CNOT gates. The explicit form of Vi in (f32|) is calculated by the 
algorism in [25| . 
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FIG. 7: The whole optical setup of simulating LHS of 4D YBE shown in Figj5] The binary numbers indicate 

location qubit basis. Each pair of —n/2 phase shifters accompanied with every beam splitter is not shown. The 

phase shifts of other PSs, from a to f, are -(n + c/>)/4, (5ir + </>)/4, n/2, {4> - 7r)/4, {4> - 37r)/4 and (n - c/>)/2, 

respectively. The angles of wave plates to form U2 and lj\, from g to 1, are n/2, (n — c/>)/8, — 7r/4, 7r/4, (57r — <fi)/8 

and 0, respectively (see ([551) ). Mirrors are placed on every corner. 



